Fractal Weyl laws for chaotic open systems 
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We present a conjecture relating the density of quantum resonances for an open chaotic system 
to the fractal dimension of the associated classical repeller. Mathematical arguments justifying this 
conjecture are discussed. Numerical evidence based on computation of resonances of systems of n 
disks on a plane are presented supporting this conjecture. The result generalizes the Weyl law for 
the density of states of a closed system to chaotic open systems. 
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The celebrated Weyl law concerning the density of 
eigenvalues of bound states is a central result in the 
spectroscopy of quantum systems Q . The Weyl formula 
states that the asymptotic level number N(k), defined 
as the number of levels with k n < k (where k — > oo) 
is given after smoothing by N(k) = #{fc„ : k n < k} = 
Vk D I (Z?/2)!(47r) £) / 2 + ..., for a quantum system bounded 
in a region R of Z?-dimensional space whose volume is V. 
For closed systems with smooth boundaries, the Weyl 
formula is well-established, and although primarily valid 
in the semi-classical limit, nevertheless can be applied 
with astonishing accuracy to very low energies extend- 
ing almost down to the ground state of integrable and 
chaotic closed systems. Generalizations of the Weyl law 
to other situations have long been sought. A notable 
example is the conjecture by Berry |2| for the density 
of states of closed systems with fractal boundaries, i.e. 
"fractal drums" . 

Open systems are characterized by resonances defined 
by complex wavevector k n = Re(fc n ) + ilm(fc„), corre- 
sponding to states with finite life times arising from es- 
cape to infinity. Open chaotic systems, which occur in a 
variety of physical situations, arc gcnerically character- 
ized by a classical phase space repeller that is fractal. In 
this letter we present a conjecture relating the density 
of resonances for an open chaotic system to the fractal 
dimension of the associated classical repeller. It can be 
stated as: 



N(k) = # {k n : lmk n > -C, Rek n < k} ~ k u 



(1) 



where du is the partial Hausdorff dimension of the re- 
peller |3|, §4.4]. This relation generalizes the Weyl law 
for the density of states of a closed system to chaotic 
open systems. 

In this letter, we will provide a heuristic argument for 
the validity of this conjecture, and present new compu- 
tations that confirm its validity. 

The repeller in a scattering problem is defined as the 
set of points in phase space which do not escape to in- 
finity at both positive or negative times. The Hausdorff 



dimension of the repeller is given by Dh = 2d# +2, where 
we did not restrict ourselves to an energy surface. For 
closed two dimensional systems, such as compact surfaces 
of constant negative curvature, we have real zeros only 
and N(k) = #{fc„ : k n < k} ~ k 2 , which is consistent 
with Q as du = 1 then everything is trapped. 

Our motivation comes from rigorous work on quantum 
resonances and in particular from the work of Sjostrand 
on geometric upper bounds on their density. The 
optimal nature of that bound was recently indicated by 
a numerical experiment |f| involving a computation of 
quantum resonances for semi-classical Schrodinger oper- 
ators with chaotic classical dynamics. 

Here we consider a different but related problem. Sup- 
pose that Z(k) is the semi-classical Selberg-Ruelle zeta 
function, with k the wave number. In some situations 
the zeros of its meromorphic continuation approximate 
semi-classically the quantum resonances of the quantized 
system - see [a, for a review of recent theoretical and 
experimental results. Because of the exact Selberg trace 
formula, this is rigorously known for surfaces of constant 
negative curvature ||. Recently that case was studied 
both theoretically and numerically in that setting. It 
was shown that a better bound is possible in small energy 
interval: # {k n : lmk n > -C, k < Rek n < fc+1} ~ k d » . 
Bounds such as these normally follow from with a 
good error estimate. Numerical results for Selberg zeta 
functions also indicate that the upper bounds are sharp. 

Here we show that {Q holds for zeta functions asso- 
ciated to configurations of hard disks in the plane. The 
result is demonstrated by numerical computations of res- 
onances of n-disks on a plane obtained from the poles 
of a semi-classical Ruelle zeta function calculated using a 
cycle expansion. We have chosen the rt-disk geometry be- 
cause it is hyperbolic, has been treated extensively with 
well-developed theoretical machinery in the form of cycle 
expansions, and has computable fractal dimensions. 

The trace formula! of Selberg, Gutzwiller, and Balian- 
Bloch (see for instance ^{j for a survey and references) 
provide one of the most elegant and useful ways of ex- 
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pressing the classical quantum correspondence. To for- 
mulate it in terms of the semi-classical zeta function 
we consider a semi-classical quantization H of a clas- 
sical Hamiltonian H: for instance H = p 2 + V(q) and 
H = -h 2 A + V{x). 

The contribution of periodic orbits (PO) to the trace 
of the resolvent is is given by 

1 , 
tr — = 

E H \p.o. ^ 

(lnZ)'(E){l + h ai {E) + h 2 a 2 {E) + •••), 
where Z(E) is the semi-classical zeta function given by 

/ °° i iu p>n +inS p (E)/h\ 

Here p's are primitive periodic orbits, T p 's their periods, 
Vp.n, the Maslov indices at n iterations, S p 's classical 
actions, and J p 's the linear Poincare maps. The formula 
is somewhat formal but can be made rigorous - see 
for a recent presentation in a generalized setting. We 
also note that, suitably modified, an exact formula for a 
class of constant curvature open chaotic systems is given 
in 0. 

For open systems, the quantum resonances are defined 
as the poles of the meromorphic continuation of the resol- 
vent (or Green's function) (E — H)^ 1 . Neglecting higher 
order terms in |J2J suggests that in the semi-classical limit 
the resonances should be approximated by the complex 
zeros of the analytic continuation of Z(E). There exists 
enough evidence now that this is the case 0. 

Since a resonance corresponding to E = Eq — iT prop- 
agates as exp(—itE/h) = exp(— itEo/h — tT / h) (hence Eq 
is interpreted as the rest energy, and T as the rate of de- 
cay), the "visible" resonances should satisfy V < Ch: if 
r h the state decays too fast to be seen. As a den- 
sity of resonance states near a given energy level, Eq, it 
is thus natural to consider NniMud) = — iT : 

\E - E \ < 5 , T < Ch} . From [1 we expect that for 
hyperbolic classical flows 

N h (E ,S)^h- d ^' 5 \ 

d(E ,S) = ^dxm{(q,p) : (4) 

\H(q,p)-E \ <5,$*(g,p) 7 4oo, t -> ±00} 

The set appearing in the definition of d(Eo,5) is the 
trapped set or the repeller of the classical flow. It is not 
clear at this point what notion of the dimension should 
be used. The upper bounds Q and numerical results in 
are given in terms of the Minkowski dimension. 

To indicate the reasons behind Q we first recall the 
standard argument for obtaining the Weyl law. If H is a 
quantum Hamiltonian with a discrete spectrum near Eq 
(which classically corresponds to the fact that H^^Eq — 



6, Eq + 6} ) is compact), we have a semi-classical trace for- 
mula H3, tvf(H) ~ hr D J J f{H{p, q))dpdq. Choosing / 
close to a characteristic function of an interval, we obtain 
the Weyl law. When the set of closed orbits has measure 
zero we have a more precise result at non-critical energies 
Eq: ixf(H - E ) ~ h- D J J f{H{p, q) - E )dpdq. 

Heuristically, the Weyl law for the number of quan- 
tum states of a confining Hamiltonian between energies 
a and b is obtained by counting the number of indepen- 
dent quantum states covering the corresponding part of 
phase space. By the uncertainty principle, a maximally 
localized state lives in a box with sides Vh (in phase 
space) and the number of such boxes needed to cover the 
region a < H < b is proportional to the volume of the 
region times h~ n where n is the number of the degrees of 
freedom. In fact, by using max-min arguments this can 
easily be made rigorous. 

For an open system the volume of the region a < H < b 
is infinite but the density of quantum resonances should 
be related to the trapped set only, that is to the set of 
points which are not carried to infinity by the classical 
flow in negative and positive directions. If that set is very 
regular in the sense of self similarity, then the number of 
maximally localized states needed to cover it is propor- 
tional to /i~ m / 2 where m is the dimension of the trapped 
set for energies between a and b. Normally we do not ex- 
pect the trapped set to be uniformly regular, especially 
when we vary the energy, but a reasonable regularity can 
be expected from its intersection with a Poincare sec- 
tion at a given energy. Then the covering argument can 
be applied on the Poincare section, and since the states 
are invariant under the flow, we obtain the missing di- 
mension. That gives the infinitesimal density at a given 
energy in terms of the dimension of the trapped set at 
that energy. 

For more quantitative variations of this argument 
(which cannot be made fully rigorous due to the non- 
self- adjointness of the problem; only upper bounds can 
be obtained following the work of Sjostrand) see || . Here 
we will present an argument based on the trace formula. 

For open systems the space in which the trace is taken 
needs to be modified. One way to do it is by con- 
jugating H with exp(G/ft) where G is a quantization 
of a Lyapunov function, (d/dt)G o $ t (g,p) ~ p{q,p) 2 , 
p{<liP) — distance to the trapped set [4|. For a suit- 
able global choice of G, the resonances are the eigenval- 
ues of the non-self-adjoint operator Hq = e~ G / h He G / h , 
and for /'s of the form exp(— az 2 — ibz), a,b > 0, 
we still have a relation between the resonances and 
tr /(Hg)- We note that the standard "Heisenberg pic- 
ture" argument and the property of G show that Hq 
is a quantization of H(q,p) — ip(q,p) 2 . Inserting this, 
somewhat formally, to the trace formula above we ob- 
tain txf(H G ) - h-° f H _ 1{Eo) exp(-p(q,p) 2 /h)dqdp ~ 
h -D h (2D-i-dimK Eg )/2 ^ waere x Eq is the trapped set in 
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H- 1 {E Q ), dimK Eo ^ 2d(E ,S) - 1. This gives the re- 
lation (@J. We stress that this argument can be made 
rigorous as far as the upper bound is concerned. 

In view of the semi-classical connection between the 
zeros of Z(E) and the resonances seen through ©, it is 
natural to consider the analogue of Q for those zeros. 
It is straightforward but a bit cumbersome (in view of 
boundaries) to rewrite the above argument in this case. 

In the case of hard disk scattering in two dimensions 
(see and references given there) the quantum Hamil- 
tonian is given by — h 2 Ao where Ac is the Dirichlct 
Laplacian. It is then natural to introduce a new vari- 
able fc, k 2 — E/ti 2 . Semiclassical asymptotics corre- 
spond to the limit k — *■ oo, and the semiclassical den- 
sity of resonances considered above should be replaced by 
N(k) = # {fc„ : lmk„ > -C, Refc„ < k} . The Hausdorff 
dimension Dh of the repeller in (0J is now independent 
of the energy level. 

For the hard disk geometry, the zeta function can be 
considered as a function of k and it takes a somewhat 
simpler form 

z(*)=«pi-x;x; — - — ^—7^1 

\ V n=l r J P I / , r s 

3=0 P 

where A p > 1 is the larger of the two eigenvalues of J p , 
and m p is the number of reflections of p. Here we set 
particle velocity v = 1 for simplicity. 

Effective ways of evaluating the analytic continuations 
of semi-classical (and dynamical) zeta functions have 
been developed by several authors. The cycle expansion 
method has proved itself to be particularly succesful. 
We used it in earlier computations performed for the pur- 
pose of comparisons with experimental data 14'j. As a 
limit case, one first considers the two-disk system. In the 
periodic orbit theory, quantum resonances form a rect- 
angular lattice in the complex fc-plane. Thus N(k) ~ k. 
This is consistent with dn = since this system is not 
chaotic. 

We also chose configurations for which the dimensions 
of the repellers were readily available three symmet- 
rically spaced disks of radii a = 1, with centers r = 6 
apart. The zeta function (|SJ with smaller j will give 
sharper quantum resonances. For example, for the sym- 
metric 3-disk system with r — 6, the quantum resonances 
for the zeta function (JSJ with j = are located in the 
area Imfc < —0.121 while the quantum resonances for the 
j = 1 zeta function are located in the area Imfc < —0.699. 
We will only consider the zeta- function with j = 0. The 
C$ v symmetry of the 3-disk system helps further factor- 
ize the zeta function and makes the periodic orbit theory 
work more efficiently |l3| . There are three irreducible 
representations, A\, A2, and E. The cycle expansion 



with PO up to period 6 gives very accurate calculation 
of the quantum resonances of the A\ and A2 represen- 
tations for the range Imfc > —0.3 and < Refc < 2000. 
For the resonances outside this range, more PO would 
be needed in the cycle expansion. For example, includ- 
ing all POs up to period 4, 5, 6, one will get respectively, 
26, 33, 39 resonances in the area 2000 < Refc < 2010 
and — 0.5 < Imfc < —0.1. The quantum resonances of A\ 
and A2 representations are distributed on lines. The E 
resonances appear to be two-dimensional, but a careful 
inspection reveals that they are also distributed on lines. 
All these lines tend to move closer to the line Imfc = — ^70 
as Refc increases, where 70 is the classical escape rate. 
This leads to the increase of the density of resonances 
0. 

FIG. 1: A2 resonances of the 3-disk system with r — 6 and 
a = l. 

The principle @ applied to the density N(k) defined 
before suggests the law 

\ogN(k)/\ogk~ D H /2 = l + d H , fc^oo. (6) 

Representative numerical results are shown in Fig^ 
for a 3-disk configuration with r = 6 and a = 1. Fig|2{a) 
shows the number of resonances N(k) vs. fc in a strip 
of width C — 0.28. Clearly the dependence of iV(fc) on 
k is superlinear. This is confirmed by the logarithmic 
plot in the bottom panel which also yields an exponent 
of 1.288, close to the calculated value of 1 + du = 1.2895. 
Fig. 2(c) shows that the density of resonances AN/ AC 
peaks at C ~ ^7o ; and then decays rapidly for large C. 
This demonstrates that the range of C used in the present 
study is more than adequate. 

Similar calculations of the quantum resonances with 
different separation r — 5, 8, 10, 12 were also carried out 
for the 3-disk system. They all confirm the validity of 
(|6"|) for the symmetric three-disk configurations. 

FIG. 2: (Top) The counting function, N(k), for width C = 
0.28 for the resonances in Fig.l. (Middle) The plot of log N(k) 
against log k. The least square approximation slope is equal 
to 1.288. (Bottom) Dependence of density of resonances 
AN/ AC on strip width C. The vertical line is ^70. 

To be valid the result of this paper requires that the ex- 
ponent of fc should be independent of the counting strip 
C. The numerical calculations indeed confirm this as 
shown in Fig[3| where only weak dependence of the ex- 
ponent on the choice of C is observed as C is increased. 

We have further checked the result by varying the disk 
separation. We also see that the agreement with the 
dimension persists when we change the disk separation 
(Fig. 3). To combine data for different r/a values in a 
single plot, we rescaled the strip width by ^70. This 
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is motivated by the fact that the density of resonances 
peaks at C ~ ^70. 

FIG. 3: Dependence of exponent on the rescaled strip width, 
2C/70, for the 3-disk system in three cases with r/a = 5, 6, 10. 
70 = 0.4703, 0.4103, 0.2802 is the corresponding classical es- 
cape rate. The solid lines are the corresponding Hausdorff 
dimensions dn=0.3189, 0.2895, 0.2330. 

The computational results presented here produce 
strong evidence for the connection between the density of 
resonant states and the fractal dimension of the repeller 
on which they concentrate. It is quite possible that dif- 
ferent dimensions might occur in lower and upper bounds 
once less symmetric examples are considered. Since zeta 
functions were used in the computation, we have also pro- 
vided evidence that the density of zeros of zeta functions 
is related to the dimension of the repeller. 

It is worth noting that while fractality arises from R 
or its boundary dR in the works of Ref. 2] , in the present 
work both R and dR are smooth and instead the classical 
phase space is fractal. 

The connections described here between quantum 
spectral properties and classical fractal properties of the 
associated repeller of open chaotic systems parallels a 
similar connection established earlier between the quan- 
tum spectral autocorrelation and the classical decay rate 
0, 0] • The spectral autocorrelation of the quantum mi- 
crowave spectra of n-disk billiards was shown to display 
the fingerprints of the classical Ruelle-P ollicott (RP) res- 
onances 7„ = 7„ ± ij' n . The leading RP resonance 70 
(70 = 0) represents the classical decay rate and is re- 
lated to the information dimension d] of the repeller by 
70 = A(l — di), where A is the average Lyapunov expo- 
nent. The higher RP resonances 7„ with n > represent 
fine structure properties of the fractal repeller comprised 
of the manifold of trapped orbits. 

In conclusion, the present work establishes a funda- 
mental connection between a quantum spectral property, 
and the fractal phase space structure of the correspond- 
ing classical dynamics. The result has important implica- 
tions for a variety of areas where it is important to know 
the number of degrees of freedom of quantum chaotic sys- 
tems, and since open quantum chaotic systems are rele- 



vant to a wide range of problems in atomic and chemical 
systems 0, lldt 1 1 5| and quantum computation. 

Work by WTL and SS supported by NSF-PHY- 
0098801. Work by MZ supported by NSF-DMS-0200732. 



* Electronic address: wtlu@neu.edu 

' Electronic address: srinivas@neu.eduj 

* Electronic address: zworski@math.berkeley.edu 

[1] H.P.Baltes and E.R.Hilf, Spectra of Finite Systems, B-I 

Wissenschaftsverlag, Mannheim (1978). 
[2] M.V. Berry in Structural Stability in Physics, eds W. 

Giittinger and H. Eikemeier, Springer, p. 51-3. 
[3] P. Gaspard, Chaos, scattering and statistical mechanics, 

Cambridge Univ. Press, 1998. 
[4] J. Sjostrand, Duke Math. J. 60 (1990), 1; M. Zworski, 

Inv. Math. 136, 353 (1999); Notices Amer. Math. Soc. 

46, 319 (1999), J. Sjostrand, in Microlocal analysis and 

spectral theory (Lucca, 1996), 377-437, NATO Adv. Sci. 

Inst. Ser. C Math. Phys. Sci., 490, Kluwer Acad. Publ., 

Dordrecht, 1997. 
[5] K. Lin, J. Comp. Phys. 176 (2002), 295; K. Lin and M. 

Zworski, Chem. Phys. Lett. 355, 201 (2002). 
[6] A.Wirzba, Physics Reports, 309, 1-116 (Feb. 1999) . 
[7] W.T. Lu, K. Pance, P. Pradhan, and S. Sridhar, Physica 

Scripta T90, 238 (2001). 
[8] S.J. Patterson and P. Perry, Duke Math. J. 106, 321 

(2000). 

[9] L. Guillope, K.Lin, and M. Zworski, to appear. 

http : / /xxx . lanl . gov/list/math . MP/new 
[10] P. Gaspard, D. Alonso, and I. Burghardt, in Advances in 

Chemical Physics XC, I. Prigogine and S.A. Rice, eds. 

John Wiley & Sons, 1995. 
[11] J. Sjostrand and M. Zworski, J. Math. Pure Appl. 81, 1 

(2002). 

[12] M. Dimassi and J. Sjostrand, Spectral asymptotics in the 
semi-classical limit, Cambridge Univ. Press, 1999. 

[13] P. Cvitanovic and B. Eckhardt, Phys. Rev. Lett. 63 
(1989), 823; P. Cvitanovic and B. Eckhardt, Nonlinearity 
6, 277 (1993); B. Eckhardt and G. Russberg, Phys. Rev. 
E 47, 1578 (1993). 

[14] K. Pance, W.T.Lu and S. Sridhar, Phys. Rev. Lett. 85, 
2737 (2000). 

[15] M. R. Haggerty, J. B. Delos, Phys. Rev. A 61, 053406 
(2000). 



-0.1 



-0.15 



-0.2 



-0.25 



-0.3 




200 



-0.1 



-0.15 



-0.2 



-0.25%* 



-0.3 




1600 



1650 



1700 

Rek 



1750 



1800 




£ o 

cu 0.1 
Q 



0.3 0.4 

Strip width C 



0.6 



0.4 



0.35 



^°Oo nO°°° 
o u o qo qu O w O 



o r .o oo00 ooo oo 



o 

3h O *^m,*** 
■ Q 

** 



°000 oUOOOOOOqOC) 



******** * *** 



o _ ; 
o 0.25 

Q. 
X 
LU 



** 



aA A ,aA a aa a ; aAa aAaa AaaA a a 



A- 



A 



A A A AA^ 



7X^ 



a a aa AaA a 



0.2 



A 



"A 



A 



0.15 



O r/a=5 
* r/a=6 
A r /a=10 



0.1 



J I I I I I I L 



1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 1.55 

strip width (2C/y ) 



